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ABSTRACT 


This  report  describes  methods  for  the  proper  design  of  first-order,  • 
recursive,  fixed-weight,  linear  filters.  Expi-essions  are  derived  and  listed 
for  commonly  used  design  parameters  such  as  noise  ratio,  transient  response, 
and  truncation  error.  The  performances  of  critically  damped  and  steady-state 
optimum  filters  are  compared.  Design  curves  are  given  that  can  bo  used  to 
select  the  weights  of  the  steady-state  optimum  filter  from  total  error 
requirements. 
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1.  Filter  Equofions 

Tlio  fii’Rt-oi’dci’  linc'ar  filters  iiiiclcr  consiclfn-atifm  are  thoae  of  the 

i'orni 

.  y  Oi  +  1  (  ■-  y  I  TV  (ii) 

y  (11  I  1)  ^  y  (n)  (1) 

y(n)  -  y  (n)  +  w,(n)  lx(n) 

y(e)  =y  (n)  -t  \V2(n)  )x(n)  -y(n)]  .  (2) 


In  state'  vec.'tor  notation,  these  equations  may  be  written 

Y(n  +  1)  =  4>(n)  Y(n) 

Y(n)  Y(n)  +  W(n){X(n)  -  HY(n)]  . 


(4) 


In  the  pi  eec'cling  equations, 


y (n  +  1) 
y(n+  1) 

is  the  predicted  system  state  made  at  time  iiT  for  time  (n  +  1)T  given n  measure 
monts.  equation 


Y(in-  1)  - 


Y{n  +  1) 


y(n) 

yfii) 


is  the  smoothed  system  state  made  at  time  nT  for  time  nT  given  n  measure¬ 
ments.  The  equation 


X(n)  |x(n)l 

l.s  the  position  measurement  corrupted  by  white  Gaussian  noise.  The  term  T 
Is  the  sampling  Intei'val,  and  n  Is  the  time  index.  TIic  equation 

T 

,1 

is  th.i  state  ti'ansition  matrix  for  a  constant  velocity  tT’ajeetory.  The  equation 

II  '  (  1  0] 


<l>(n)  - 


I 


.  {  ' 


is  the  observation  matrix.  The  weighting  matrix  is 

wj(n) 

■  :  '  W(n)  =  W2(n) 

In  the  application  considei'cd  for  this  filter,  y  and  y  are  estimated  posi¬ 
tion  and  velocity,  and  x  is  measured  xjosition  along  one  coordinate  axis/  Equa¬ 
tion  (1)  or  (3)  is  called  the  prediction  equation,  and  Equation  (2)  or  (A)  is  called 
the  smoothing  or  correction  equation.  The  objective  of  this  paper  is  the 
selectj^on  of  a  constant  weighting  matrix  \V  to  get  the  best  estimate  of  the  predicted 
state  Y.  In  m(^t  of  the  literature  dealing  with  a  constant  weighting  matrix, 

Wj  is  called  g  or  a,  and  TW2  is  called  h  or  0.  The  g-h  notation  will  predominate 
in  this  paper. 

The  overall  filter^ performance  may  be  described  by  tliree  factors.  They 

•  are: ,  ,  .s/' 

a)  Variance  reduction  ratio 

b)  Transient  response 

c)  Truncation  error." 


A  fourth  and  sometimes  dominant  factor  is  computional  complexity.  In  tracking 
many  targets  by  use  of  a,  time  shared  computer  and  radar,  this  factor  limits 
the  filter: order  to  first  or  second,  necessitates  the  use  of  constant  weights, 
and  makes  smoothing  in  one  coordinate  at  a  time  highly  desirable. 


2.  Voriance  Reduction  Ratio 

The  variance  reduction  ratio  is  the  ratio  of  the  rms  noise  output 
from  the  filter  to  the  rms  noise  input.  It  is  sometimes  more  appropriately 
called  the  noise  amplification  factor,  since  it  may  necessarily  exceed  unity  if 
the  filter  requirement  is  for  good  dynamic  response.  In  general,  low  noise 
output  requires  small  weights,  and  good  transient  performance  requires  large 
weights.  A  mathematical  expression  can  be  derived  for  the  first  order  filter 
in  which  Wi=  g  and  W2  =  h/T  are  constants.  If  equations  (1)  and  (2)  which  I 

occur  at  time  nT  are  combined  with  the  same  equations  for  time  (n-l)T.  all  v 

quantities  except  predicted  positions  and  measurements  may  be  eliminated  to  I 

yield  '  •; 

y(n+  1)  =  fg+ h]  x(n)  -  gxfn  -  1)  +  [2  -  g  -  h]  y(n)  j 

+  [g  -  1]  y(n  -  1)  .  . 

'  '  h 

..  ■ '  f 

2 

I 

jf 

s 


r 

I  * 


Equation  (5)  gives  predicted  position  in  the  standard  feedlDaclc  filter  form 

y(n  +  1)  =  ^  a.  x(n  -  i)  -  X 

i=0  '  3=0  J 

In  Equation  (5) ,  the  output  is  the  weighted  sum  of  the  last  two  inputs  and 
outputs. 

Since  the  filter  is  linear,  the  Superposition  Theorem  may  be  invoked 
in  calculating  the  noise  ratio  which  is  independent  of  any  trajectoiy  dynamical 
errors.  If  the  input  is  sampled  white  noise  with  variance  =  E  [x"] ,  in 

which  Ef  ]  denotes  "expected  value  of,"  then  the  output  noise  Variance  is 

Efv^(n+  1)].  Squaring  Equation  (!5)  and  taking  the  expected  value  of 

y  ‘  , 

the  result  gives 

(g+ h)^/^/+ (2-g-h)2cr-2+  (g-l)^o*^ 

_  2g  (2-g-h)  E|x{n- 1)  v(n))  +  2  (2-g-h)  (g- 1)  Efy(n-l)yfn)]  . 

.i.:  (7) 

In  arriving  at  Equation  (7),  it  was  assumed  that  the  input  and  output  noise 
distribhtions  arc  stationary  so  that 

‘  E|v2(n^  1)1  =  EFy2(n)]  =  Efy2(n  -  1)]  = 

and 

E I  x‘ (n )  ]  =  E  ( x^  (n- 1 )  1  =  . 

Also,  the  input  noise  samples  arc  independent  so  that 

..  Efx(n)  x(n  -  1)  ]  =  0 

and  pi’cdictidns  arc  independent  of  future  measurements  so  that 

Elx(n)  y(n)l  Elx(n)  y(n  -  1)J  =  E|x(n  -  1)  y(n  -  1)1  =  0. 

The  remaining  two  terms  in  Equation  (7)  involve  correlated  quantities  but  can  be 
evaluated  from  Equation  (.5).  Writing  Equation  (5)  fory(n),  multiplyingby 
X  (n  -  1) ,  and  taking  the  expected  value  of  the  result  gives 

E|xfn  -  1)  y(n)]  •=  (g  +  h)  ^ 

•I 


injphich  use  has  been  made  of  the  same  stationarity  and  independence  assumptions 
a^before.  Again,  writing  Equation  (5)  fory(n),  multiplying  by  y  (n-1) ,  and  . 
tafang  the  expected  value  gives 


E[y(n)  y(n-l)]  = 


(2.-g-h)  CT-/  -  g(g+  h)  (j  ^ 

y 


(9) 


after  making  the  same  assumptions  as  before  and  making  use  of  Equation  (8). 
Fmally,  substituting  Equations  (8)  and  (9)  into  Equation  (7)  gives 


■) 

y  _  2g-  gh  +  2h 

(j^2  g(4-2g-h) 


(10) 


Equation  (10)  gives  the.  i*atio  of  the  variance  of  output  noise  in  predicted  posi¬ 
tion  to  that  of  measurement  noise  for  constant  filter  weights.  This  is  one 
equation  applicable  to  making  a  proper  choice  of  g  and  h.  Similar  equations 
ca^bc  derived  for  smoothed  position  and  predicted  velocity.  They  arc: 


^y  2g^  +  2h  -  -Sgh 
“  g(4-2g-h) 


(11) 


^y  _  1  2h^ 

^  g(4  -  2g  -  h 


(12) 


Equations  (10),  (11),  and  (12)  can  also  be  derived  from  the  system 
unit  impulse  response  by  using 

oo 

^2=  E  '  (19) 

X  0 


in  which  h(n)  is  the  inverse  transform  of  H(Z)  =  Y  (Z)/X  (Z) ,  and  Y  (Z)  is 
obtained  by  taking  the  transform  of  Equation  (.5)  for  predicted  position. 
Equation  (13)  can  be  evaluated  from  the  integral  for  the  sum  of  a  squared 
sample  sequence  [1]. 


Eh2(n)=^  ^H{7.)  H(z-')  7-‘‘ dz  ,  (14) 

0 

in  which  the  path  of  integration  is  the  unit  circle  in  the  Z^planc,  The  mean 
square  noise  ratio  then  becomes  27rj  times  the  sum  of  the  residues  of 
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Equation  (14)  at  the  interior  poles  of  the  integi'and.  However,  it  is  felt  that 
the  first  method  of  derivation  more  clearly  points  out  the  assumptions  that  have 
to  be  m.ade  in  deriving  Equations  (10),  (11),  and  (12).  The  same  assumptions 
have  to  be  made  in  deri\'ing  Equation  (13). 

The  assumptions  of  stationarity  should  be  kept  in  mind  when  applying  equa¬ 
tions  such  as  Equation  ( 10) .  The  input  Stationarity  assumption  merely  implies  that 
the  root  mean  square  measurement  noise  does  not  change  significantly  during 
one  sample  period.  This  is  entirely  reasonable  from  the  standpoint  of  the 
trajectory  itself;  that  is,  measiu'ed  coordinates  do  not  change  fast  enough 
between  samples  to  affect  the  expeeted  noise,  which  is  a  function  of  position. 
The  change  in  measurement  noise  due  to  target  scintillation  is  another  matter. 
In  practice,  the  target's  effective  radar  cross  section  can  vary  with  time  such 
that  significant  de\iations  i.i  the  svstem  signal-to-noise  ratio  oc  -ur  at  rates 
close  to  normally  used  tracking  rates.  However,  there  is  evider  ce  that,  for 
periodic  variations  in  signal  to  noise  in  the  steady  state.  Equation  (10)  still 
gives  the  ratio  of  total  mean  square  noise  output  to  total  mean  sq' are  noise 
input  [2].  Even  though  the  input  noise  is  stationary,  the  output  noise  need  not 
be.  There  is  a  transient  peidod  following  track  initiation  in  the  noise  ratio 
itself  about  which  not  enough  is  known  at  this  time.  One  reference  indicates 
that  Equation  (10)  gives  an  entirely  erroneous  result  during  this  transition 
period  1 3] .  Further  work  is  planned  in  this  area. 


3.  Transient  Response 

The  specification  of  system  transient  response  has  always  been 
somewhat  of  a  problem.  The  classical  method  is  the  use  of  time  constant,  or 
time  for  the  output  to  decay  to  l/e  of  some  initially  stored  value.  Other  criteria 
such  as  10  to  90  percent  rise  time  have  been  used.  These  ideas  are  still 
applicable  as  definitions  of  transient  response  in  sampled  data  systems,  but  a 
more  elegant  approach  is  to  calculate  the  sum-squared  system  error  output 
ca,used  by  some  input  for  which  the  output  should  eventually  converge.  As  was 
indicated  in  Equation  (14),  if  f(nt)  is  a  convergent  series  of  data  samples 
whose  ultimate  value  is  zero,  then  the  sum  of  the  squares  of  all  the  samples 
from  zero  to  infinity  is  given  by 

~  ^  F(z)  F(z”’)  z”'dz  .  (15) 

ZTT] 

In  Equation  (15),  F(z)  is  the  z  transform  of  f(nt)  defined  by 


F(z)=J^.  f(nT)z"" 
n=0 


(Ifi) 


o 


and  the  path  of  integration  is  the  unit  circle. 

The  transform  of  Equation  (5)  is 

O,  .  g  +  h  -  g  z-  ^ , 

1-  (2-g*h)z  (g-Dzr" 

Since  the  first-order  filter  will  converge  on  a  position  ramp  input  (velocity 
step)  without  error,  a  reasonable  input  to  be  used  as  a  basis  for  transient 
response  calculation  is 

x(nT)  =  nT  ,  ( 


whose  transform  is 


X(z)  = 


Tz~’ 
(l-z-’)  = 


The  desired  prediction  is 


(19) 


Y(z)  =  zX(z)  =  •  <20) 

The  transform  of  the  error  is 

e-(z)  =Y(z)  -Y(z)  =  i_(2-g-h)z-' -  (g-l)z-“ 
after  substituting  Equations  (15),  (17),  and  (18).  Evaluating  the  integral 


by  summing  the  residues  at  the  two  poles  inside  the  unit  cii’cle,  gives 


T^(2-g) 

gh(4-2g-h) 


(23) 


for  the  sum-squared  error  in  predicted  position  of  the  first  order  filter  with 
a  ramp  input.  Similarly,  it  can  be  shown  that  for  smoothed  position  and 
velocity. 
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4.  Truncation  Error 


gh(4-2g-h) 

Vi_g^(2-g)4  2h(l-g 
gh  (4  -  2g  -  h) 


It  was  stated  in  the  previous  section  that  a  first-order  filter  will 
converge  on  a  constant  velocity  noise  free  input  trajectory  with  zero  error.  It 
is  also  true  that  a  constant  acceleration  input  will  be  estimated  with  a  constant 
error  in  the  steady  state,  and  trajectories  with  higher  derivatives  will  cause 
errors  that  grow  as  long  as  the  higher  derivatives  do  not  change  sign. 
Fortunately,  real  target  trajectories  cannot  possess  monotonic  higher 
derivatives  along  any  given  coordinate  axis  for  verv  and  the  position  errors 
that  result  can  be  kept  within  reasonable  limits.  The  steady-state  error  that  is 
caused  by  higher  derivatives  than  the  first  is  called  truncatibrl  error  for  the 
first-order  filter. 


The  mathematical  approach  to  estimating  truncation  error  is  to  assume 
that  that  due  to  acceleration  is  predominant.  This  is  equivalent  to  assuming  that 
acceleration  does  not  change  significantly  in  one  sample  period.  This  is  a 
good  assumption  for  trajectories  for  which  a  first-order  filter  is  useful.  In 
practice,  it  is  found  that  tracking  filters  higher  than  second  order  are  rarely 
necessary  and  a  first-order  filter  is  frequently  adequate.  It  is  also  true  that 
the  noise  amplification  ratio  tends  to  increase  with  the  order  of  the  filter; 
that  is,  a  low-order  filter  will  follow  a  low-order  curve  with  less  noise  output 
than  a  higher  order  filter.  The  order  of  the  filter,  then,  like  the  choice  of 
smoothing  constants  is  a  compromise  between  noise  reduction  and  dynamic 
response. 

The  fact  that  truncation  error  for  the  first-order  filter  can  be  closely 
approximated  for  several  representative  trajectories  by  the  response  to  the 
second  derivative  input  has  been  demonstrated  experimentally  [4] .  The 
anaMical  investigation  can  proceed  as  follows :  Let  the  input  be 


x(nT)  =  A./2(nT)^  ,  (26) 

which  corresponds  to  a  parabolic  trajectorv'  in  one  coordinate  with  acceleration 
A.  The  2  transform  is 


(27) 
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5.  Wftight  Selection 

The  weights  g  and  h  nuist  be  selected  te  give  some  suitable 
eoiii[)romise  between  dynamic  rosponpe  and  noise  ratio  which  meets  the  overall 
system  error  specification.  Sometimes  the  sampling  period  T  can  btj  varied 
within  limits  to  aid  in  meeting  the  requlrCmt  ats.  Two  basic  relationahips 
between  g  and  h  have  been  used  to  ease  the  task  of  choosing  the  witights  from 
all  possible  g-h  ixtirs. 


a.  yhe  Critically  Dampj^fl  Filter 

As  can  be  seen  from  Equation  |17),  the  first-order  ttlt  -r  has  a 
soeond-degree  polynomial  in  its  transfer  function.  Ithaa  beeii  recogni/jcd  for 
many  years  that  the  critically  clamped  condition  in  a  syatem  charactnrized  by 
such  a  funcitlon  has  reasonably  fast  and  Wei l-behavad  transient  response.  The 
poles  of  Equation  (17)  are 


Zi.j  =  4  l2-g-h  iNr(2r"g-h)?+ 4(g^i)  ]  ,  '  (35) 

and  critical  damping  occurs  when  .  , 

(2-g-h)2+ T(g-l)  0  (36) 

from  which  ... 

h  =  2-g  ±  2A/TIg  .  (37) 

AlaOt  by  applying  Uouth’a  method  to  the  denominutor  polynomihal  of  Eqwtion 
(17)  after  using  th(!  bilinear  ti’ansiformation  [.5], 


S+J_ 

s"'-*  1  •  ■ 


The  lirst-brden’  filter  is  stable  If 


h  >  0 
g  >  <■ 

!-2g  -h  >  o  .  .  (39) 


« 


The  conditions  in  Equation  (39)  indicate  that  g  and  h  lie  within  the  triangle  shown 
in  Figure  1.  The  locus  of  Equation  (37)  is  also  shown  in  Figure  1  showing. the 
combinations  of  g  and  h  normally  used  for  the  critically  damped  filter.  Actually, 
the  useful  part  of  the  curve  lies  between  h  =  o  and  h  =  1,  since  for  larger 
values  both  noise  ratio  and  time  response  increase.  For  example,  if  the  time 
response  is  defined  to  be  the  l/c  time  constant  in  position  error  which  results 
from  a  velocity  step  input,  the  inverse  transform  of  Equation  (21)  with  equal 
denominator  roots  gives  for  the  time  constant  the  value  of  n  for  which 

log^  i  (n  +  1)  (1  -  g)  n/2l  =  -  1  .  (40) 

If  Equation  (40)  is  plotted  against  the  noise  ratio  [Equation  (10)  ] ,  Figure  2  is 
obtained.  The  lower  portion  of  the  curve  is  the  useful  range  and  is  obtained 
using  the  mi,ius  sign  in  Equation  (37).  The  upper  curve,  obtained  using  the  plus 
sigh,  gives  more  noise  for  the  same  time  response. 


'  b.  The  Optimum  First-Order  Filter 

There  are  two  weights  to  be  determined  in  the  first-order 
filter,  and  specif5dng  one  of  them  is  equivalent  to  fixing  one  performance 
criterion,  such  as  noise  ratio.  Then  a  relationship  exists  between  g  and  h  which 
minimizes  the  transient  error  for  this  specified  noise  ratio.  This  means  that 
the  function 

6-2+ X  22..  (41) 

■  '  y  y 

in  which  X  is  a  constant  and  6^  and  2^  are  given  by  Equations  (10)  and  (23) 

must  be  differentiated  with  respect  to  g  and  h,  set  equal  to  zero,  and  h  solved 
for  in  terms  of  g.  The  result  is 

h=g2/(2-g)  .  (42) 

Further  details  are  given  by  Benedict  and  Bordner  [6].  It  might  also  be  pointed 
out  that  Benedict  and  Bordner  showed  that  not  only  does  Equation  (42)  minimize 
Equation  (41)  but  that  the  filter  given  by  Equation  (5)  is  the  optimum  linear 
filter  for  tracking  a  noisy  ramp  among  all  the  filters  of  the  form  in  Equation  (6). 
That  is,  no  other  combination  of  past  measurement  and  predictions  will  do  any 
better  than  the  last  two. 
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V-  * 


Ktij>iiui>  3lt)tnl,v“Siuic!  1' iitei' 


II'  the  Kalman  filtering  technique  1 7  j  is  applied  to  the  estimation 
of  a  constant  velocity  trajectory  given  noisy  position  measuiunnonts,  Equations 
(.1)  and  {4)  give  the  predicted  and  corrected  states.  The  weights  are  competed 
for  each  measurement  from 


W(n)  =  P(n)  II'^  (HPH^+  R]-’ 
P  (n)  'l>P(n  -  1)  +  Q 


P  (n-1)  =  (I  -  W(n-l)  HI  P(n-l)  , 


(43) 


in  which  P  is  the  covariance  error  matrix  in  predicted  state  and  P  is  the  corrected 
covariance  error  matrix.  Other  matrixes  were  defined  previously  except  for  Q, 
which  will  be  taken  to  be 


Q.- 


0 


0  (T 


M 


(44) 


In  defining  Q  in  this  manner,  all  deviations  in  the  trajeetoi^  caused  by  deriva¬ 
tives  higher  than  the  first  are  treated  as  white  noise  with  variance*  and  this 
maneuver  noise  is  assumed  Iq  affect  only  the  unceitainty  in  predicted  velocitj' 
during  the  suinplc  period  over  which  it  occurs. 


Expiinding  the  matrixes  in  Equation  (4.T)  leads  to  the  following  equations 
for  the  elements  of  VV<n)  and  P('n): 


Wi(n)  P]i(n)/1  Pjl(ii)  +  K] 

W2(n)  =  Pi2(n)/(P„(r')  +  H] 


(45) 


1,(11)  =  W,(n-1)  U+  2T  Wsdi-l)  K  +  P22(n~l)  -  T'^  W2(n-1)  P|2(n-1) 


]n(n)  "  W2(n-1)  R+  T  Pjopi-l)  -  T  \V2(n-I)  P,2(ri-1) 


P32(n)  =  P22(n-1)  -  Wjtn-l)  p'ljdi-])  +  M  , 


(4(1 ) 


in  whi(h 


It  -  <T  ^  is  the  niensuronicnt  variance  and 

X 

M  =  is  the  inarieuver  variance  .  (47) 

M 


If  K  and  M  are  constant,  Vi^  (n)  and  P(n)  become  constant  matrixes  as  the 
time  index  n  increases.  In  the  steady  state,  the  time  Indexes  In  Eriuation  (4fi) 
may  be  dropped  and  the  constant  matrix  cleiments  solved  for.  If  this  is  done, 
there  is  obtained 


from  which 


or  in  g-h  notation, 


P',1  -  T  9,0  =  W,K  +  T  WjR 


W2  "P,2  «  M  , 


TW2  -  W,V  (2-W,) 


M 


VV, 


r“  (2-w,)'5  (i-vv,) 


(4«) 


(49) 

(50) 


,  h  =  gV  (2-g) 


(51) 


_ s'* 

(2~g)^  (l-g) 


.52) 


Equation  (51)  is  the  same  as  Equation  (42),  which  was  derived  by  calculus  of 
variations.  Thus,  this  steady-state  filter  is  the  same  as  Uenedict  and 
Bordner's  optimum  fixed  wfeight  filter.  Ildwever,  the  Kalman  filter  fixes  the 
weights  through  Equation  (.52)  which  gives  g  In  terms  of  the  sampling  period* 
measurement  noise,  and  expected  maneuver  noise.  In  the  optimum  fixed  weight 
fiitcr,  the  final  choice  of  weights  ts  still  a  compromise  between  Equations  (10) 
and  (23).  The  steady-state  Kalman  makes  the  choice  automatically  if  values 
can  be  selected  for  T,  M,  and  R. 


6.  Filter  Ccraporison  and  Design  Curves 

One  objective  of  this  paper  is  to  compare  the  critically  damped 
Alter  with  the  steady-state  Kalman  filter.  Some  basis  for  compai’laon  must  be 
chosen  and  that  basis  will  be  one  of  relative  transient  and  truncation  errors  for 
equal  noise  ratios.  Equations  (37)  and  (42)  have  been  plotted  In  Figure  3 
shpwing  h  in  terms  of  g  for  the  two  filters.  In  Figure  4,  Equations  (10),  (23) , 
and  (31)  for  noise  ratio,  transient  error,  and  Irimca'ion  error  are  shown 
(normalized)  for  the  critically  damped  filter  as  functions  of  g.  The  same 
quantities  are  shown  for  the  ateady-slate  filter  in  Figure  fi.  In  Figure  C,  the 
sum-squared  transient  errors  and  truncation  errors  for  both  filters  are  shown 


di 

t 

r- 

V 


runi;iion,s  of  niean-Hquare  nolHO  t  atio.  Figuro  fi  show  s  tltat  for  a  given  oiifpnf 
noino  irii'Ao  the  steaily-atute  filter  has  a  tranHlont  error  about  Hfi  percent  as 
great  uh  that  of  the  critically  clampccl  filter  and  a  truncation  error  aliout  72 
perceoi  iiM  great i 


In  Figures  7  (a)  and  (b) ,  g  and  h  are  shown  as  functions  of  the  parameter 
T^M/R  for  the  steady-state  filter.  Those  figures  cover  a  wide  range  because  M 
is  normally  a  function  of  so  that  T^M/R  is  a  function  of  T'*.  In  general, 

Figure  7(a)  is  useful  for  tracking  Intervals  on  the  order  of  a  second  and  Figure 
7(b)  for  Intervals  around  0.  1  second. 


Figures  7  (a)  and  (b)  are  valid  regardless  of  how  M  is  chosen.  In 
practi<!c,  T  (or  its  allowable  range)  is  known  and  a  good  estimate  can  be  made 
of  R,  the  measurement  variance,  A  reasonable  way  to  estimate  M  is  to  assume 
some  maximum  acceleration,  for  the  trajectoiT  so  that  the  greatest 

veloeity  change  u  one  sa' iple  period  is  TA  .  Let  this  quantity  be  related 
.  _  nicix 

to  a  by 
n 


=  TA 


max 


(53) 


so  that 


and 


or 


2 

max 


(5i) 


T<  A^ 


max/  X 


(55) 


IT 

X 


A 


max 


Bcr 


(56) 


If  the  maneuver  were  truly  random  from  sample  to  sample,  B  would  be  set 

equal  to  3  so  that  TA  W'ould  be  the  Ikr  value  of  the  maneuver  noise.  But  a 
max 

velocity  maneuver  is  not  really  white  noise.  If  the  target  accelerates  at  all, 
there  will  be  a  high  degree  of  correlation  between  samples  in  the  maneuver. 
There  are  techniques  available  for  mathematically  modeling  the  maneuver  to 
aceount  for  this  correlation,  but  this  introduces  more  states  to  be  estimated. 
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Like  the  question  of  track  initiation,  this  1h  u  Hubjeet  best  loft  for 
future  investigation.  Some  aid  in  selecting  I)  can  bo  obtained  by  plotting  curves 
of  truncation  error  and  noise  ratio  as  in  Figure  H.  The  normalized  I'itT  noise 
error  is  .‘id  from  Equation  (10)  which  is  plotted  as  a  function  of  g  using  Equation 
(43)  for  h.  The  normalized  maximum  truncation  eri’or  can  be  obtained  by 
equating  Equation  (!jG)  to  the  squama  root  of  Equation  (132)  ; 


T^  A  2 

_ max  _ 

.110-^  "  "  (i-g)  -  g 


Recognizing  that  the  truncation  error  is 


A 

max 


A 

h 


mid 


h  =  gV  (2-g)  , 


then  the  normalized  maximum  truncation  error  is 


A 

max  B 
0  "  \l  l-g' 


(b7) 


(58) 


(59) 


((10) 


Equation  (61)  is  plotted  as  a  function  of  the  parameter  B  in  Figure  8  and  shows 

that  the  filter  weights  the  truncation  and  noise  errors  equally  for  B  between  1 

and  2.  In  Figure  9,  B  is  plotted  against  the  worst  expcctrjd  error, 

3(5+  A  /cr  ,  with  T^'  A  Ar  as  a  parameter.  From  these  curves,  an 
max/  X  max/  x 

estimate  can  be  made  of  the  v?  i  -e  of  B  which  ininirnizca  the  maximum  error. 

In  Figure  10,  the  total  error  is  plotted  as  a  function  of  T*  A  hr  with  B  as  a 

max/  X 

parameter.  Elgure  9  is  useful  when  the  absolute  minimum  total  error  must  bo 
obtained  In  order  to  maintain  track.  Figure  10  is  more  useful  against  maneuver¬ 
ing  targets  when  the  least  obtainable  error  Is  not  required. 


Also  in  Figure  10  is  shown  the  total  error  for  the  case  g  =  1,  h  =  1,  or  no 
smoothing,  given  by 


y  (n+  1)  =  2x(n)  -  x(n  -  1) 


for  which 


(5  =  'v/IT and  A  «  T^A  . 


(61) 


(62) 
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It  can  be  dodueoci  from  Figure  9  tliat  when  truncation  error  can  be 

expected  to  exist  throughout  the  trajectory,  such  as  when  tracking  a  ballistic 

missile,  the  choice  of  D  =  1  or  leas  is  optimum.  However,  for  small  values  of 

T*'  A  h  ,  B  loss  than  one  does  not  drastically  change  the  total  error  ii  om 
ma::  /  x 

what  It  is  with  II-  1,  It  is  further  seen  from  Figure  10  that  for  ratios  of 

T“  A  h  greater  than  4,  the  reduction  in  error  over  that  of  the  no-feedbaok 
max/  X 

filter  is  eompleteiy  negligible  and  a  no-fecdback  filter  should  be  used. 

Against  a  maneuvering  target,  the  lai’gost  value  of  B  should  be  used 

Such  that  the  worst  expected  error,  A  r  3Ar  ,  does  not  exceed  the  confines 

max  X 

of  the  radar  measurement  volume.  This  will  insure  that  track  is  not  lost 
during  a  maximum  maneuver  while  giving  the  least  noise  during  no  maneuver. 

Incidentally,  Figure  10  docs  not  indicate  that  values  of  B  to  the  right  of 
the  no  smoothing  curve  arc  useless.  Aa  B  decreases,  the  cuiwes  move  across 
the  page  to  the  left  then  back  toward  the  no  smoothing  lind  as  B  gets  very  small. 
No  smoothing  corresponds  to  13  =  0. 

As  an  example  of  the  use  of  the  curves,  suppose  that  A  h  is 

estimated  to  be  0.  1  and  the  normalized  coordinate  extremity  of  the  radar 
volume  is  O.O.  From  Figure  9,  the  optimum  value  of  B  Is  about  0.6  and  the 
worst  error  is  .'t.  0,  well  within  the  radar  volume.  The  noise  ratio  from 
Figure  H  is  2.  15.  If  the  radar  volume  dimension  C.  0  is  used  for  A  h  +  .‘(6, 

tlien,  from  Figure  10,  B  4  and  the  noise  froiii  Figure  8  is  1.22.  Thus  the 
no-maneuver  error  is  about  half  what  it  la  with  the  least  maximum  error  condi¬ 
tion,  l)ut  track  is  maintained  in  cither  case. 


7.  Conclusion 

Tlie  mathematical  techniques  useful  in  the  analysis  of  first-order 
digital  filters  have  been  review'ed  and  demonstrated.  The  impoitant  performance 
measures  (noise  I’atlo,  transient  cri’or,  and  truncation  error)  have  been 
deiived  and  listed  for  predicted  position,  smoothed  position,  and  veldelty.  The 
assumptions  made  in  those  derivations  have  been  noted  and  their  efl'ecta  dis- 
cuHsed  whore  possible.  The  steady-state  optimum  filter  has  been  compared  to 
the  critically  damped  filter  and  shown  to  have  less  truncation  and  transient 
error  for  u  given  noise  ratio.  Design  curves  have  been  provided  that  can  be 
used  to  selcet  the  smoothing  weights  if  the  measurement  noise  and  trajeetoi’y 
dynamics  can  be  estimated.  This  filter,  however,  in  which  tho  weights  are 
constant,  Is  not  presently  reeommiended  for  track  initiation. 


( 
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For  ballistic  missile  tracking,  the  truncation  and  noise  errors  should  be 
weighted  about  equally  with  the  parameter  B  =  1.  A  larger  value  of  B  would  be 
better  against  a  maneuvering  target,  provided  the  total  error  does  not  exceed 
the  radar  measurement  volume. 


FIGURE  2.  SMOOTHING  VERSUS  TIME  CONSTANT  FOR  CRITICALLY  DAMPED  g-h  FILTER 


TRANSIENT  AND  TRUNCATION  ERROR 


10,0001“  100 


S  10O 


1.0  •-  0.01 


FIGUUE  5.  NOISK  IIATIO  TlUNSIENl'  EimOll  AND  TRUNCATION  ERROR 
FOR  THE  OPTIMUM  FILTER 


WEIGHTS  g  AND  h 


a 

FlGUllE  7.  K  AND  h  AS  A  FUNCTION  OF  T 


FIGURE  8.  MAXIMUM  NOISE  AND  TRUNCATION  ERRORS  WITH  B 
AS  A  PARAMETER 
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